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The current status of carlomat, a program for automatic computation of the lowest order cross sections of 
multiparticle reactions is described, the results of comparisons with other multipurpose Monte Carlo programs 
are shown and some new results on e + e~ — > bbbbuddu are presented. 



1. MOTIVATION 

Many interesting aspects of the Standard 
Model (SM) and models beyond it can be stud- 
ied through investigation of reactions involving 
a few heavy particles at a time. Owing to the 
high energy and luminosity such reactions will be 
copiously observed at particle colliders such as 
the Large Hadron Collider (LHC), or the Inter- 
national Linear Collider (ILC) PQ. As the heavy 
particles are usually unstable, they almost im- 
mediately decay leading to reactions with several 
particles in the final state. Already in the low- 
est order of SM matrix elements of such multi- 
particle reactions receive contributions typically 
from many thousands of the Feynman diagrams, 
most of which constitute background to the "sig- 
nal diagrams" representing the interesting sub- 
processes of production and decay of those heavy 
particles. Because of the large number of Feyn- 
man diagrams involved, reliable SM predictions 
for such reactions can be obtained only through 
a fully automated calculational process. 

To be more specific, let us consider, e.g. a re- 
action of associated production of the top quark 
pair and Higgs boson at the ILC 

e + e~ -> ttH. (1) 
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Because its cross section is by far dominated by 
the Higgsstrahlung off the top quark line, reac- 
tion ([T]) can be used to measure the top-Higgs 
Yukawa coupling [2. . As the top and antitop de- 
cay, even before they hadronize, predominantly 
into bW + and bW~, respectively, and the Higgs 
boson, dependent on its mass nin, decays mostly 
cither into a 66-quark pair or an clectroweak 
(EW) gauge boson pair and the EW bosons sub- 
sequently decay, each into a fermion-antifermion 
pair, reaction (QJ will lead to reactions with ei- 
ther 8 or 10 fermions in the final state. If 
ma < 140 GeV, which is favoured by the direct 
searches at LEP [5] and theoretical constrains in 
the framework of SM [3], then the Higgs boson 
would decay mostly into a 66-quark pair resulting 
in reactions of the form 

e+e- -» bbbbfj[hf> 7 (2) 

where fi, = v e ,v^,v T ,u,c and f[,f 2 = 
e~, /i~ , t~ , d, s are the decay products of the W- 
bosons coming from decays of the t- and i-quark. 
Thus, in this case reaction (JTJ) can be detected 
in any of the following channels: the leptonic, 
semileptonic and hadronic channels, as repre- 
sented by the reactions 

e + e~ — > bbbbr^ v T yT 9^, (3) 
e + e~ — > bbbbudfi" D^, (4) 
e + e~ — > bbbbusdc, (5) 

respectively. Taking into account both the EW 
and quantum chromodynamics (QCD) lowest or- 
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der contributions in the unitary gauge, with the 
neglect of the Yukawa couplings of the fermions 
lighter than c quark and r lepton, there are 
21214, 26 816 and 39 342 Feynman diagrams of 
reactions and ([5]), respectively. If both 

W + and W~ decay into quarks of the same fam- 
ily we obtain, e.g. the reaction 

e + e~ — > bbbbuddu (6) 

which, neglecting the light fermion masses, re- 
ceives contributions from 185 074 Feynman di- 
agrams. Most of the diagrams comprise back- 
ground to the 20 signal diagrams representing res- 
onant production and decay of the top quark pair 
and Higgs boson. For illustration, the represen- 
tative signal diagrams of reaction ([6]) are shown 
in Fig. [TJ 

There exists several multipurpose Monte Carlo 
(MC) generators such as HELAC/PHEGAS [5], 
AMEGIC++/Sherpa [6], 'Mega/Whizard [7], 
MadGraph/MadEvent [8], ALPGEN [9], or 
CompHEP/CalcHEP [10]. However, one may en- 
counter problems while trying to obtain reliable 
SM predictions for multiparticle reactions as ([3])- 
([5]), not to mention reaction ©, with publicly 
available versions of the generators. 

In this lecture, the current status of carlomat, 
a new program for automatic computation of the 
lowest order cross sections of multiparticle reac- 
tions is described, the results of comparisons with 
other multipurpose MC programs are shown and 
some new results on reaction ([6]) are presented. 

2. A PROGRAM 

carlomat is a program written in Fortran 
90/95. It generates the matrix element for a user 
specified process and phase space parametriza- 
tions, which are later used for the multichannel 
MC integration of the lowest order cross sections 
and event generation. The program takes into ac- 
count both the EW and QCD lowest order con- 
tributions. Particle masses are not neglected in 
the program. The number of external particles is 
limited to 12 and only the SM is currently imple- 
mented in the program. 

carlomat works according to the following 
scheme. User specifies the process he wants to 



have calculated. Then topologies for a given 
number of external particles are generated and 
checked against Feynman rules which have been 
coded in the program. In this process, helicity 
amplitudes, the colour matrix and phase space 
parametrizations are generated. Finally, they are 
copied to another directory where the numerical 
program can be executed. 

2.1. Generation of topologies 

Let us consider models with triple and quartic 
couplings. The process of generation of topologies 
starts with 1 topology of a 3 particle process that 
is depicted in Fig. O 




Figure 2. The single topology of a 3 particle pro- 
cess. 



The 4 topologies of a 4 particle process which 
are depicted in Fig. [3] are obtained by attaching 
line No. 4 to each line and to the vertex of the 
graph in Fig. [2] Then the 25 topologies of a 5 
particle process are obtained by attaching line No. 
5 to each line, including the internal ones, and to 
each triple vertex of the graphs in Fig. [3J 

The number of topologies grows dramatically 
with the number of external particles. For a 



No. of particles No. of topologies 



6 220 

7 2 485 

8 34300 

9 559405 

10 10 525 900 

11 224449 225 



process with n external particles, topologies for 
ri — 1 particles are generated recursively and 
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(a) (b) (c) 

Figure 1. Representative signal Feynman diagrams of reaction ([6]) in the unitary gauge. The remaining 
diagrams are obtained by all possible permutations of the two b and two b lines. The Higgs boson coupling 
to electrons has been neglected. 




Figure 3. The topologies of a 4 particle process. 



then, while adding the n-th particle, the program 
checks whether a topology results in a Feynman 
diagram or not. Topologies can be generated and 
stored on a disk prior to the program execution. 

2.2. Feynman diagrams 

Actual external particles are assigned to lines 
l,2,3,...,nina strict order. Each topology is di- 
vided into two parts which are separately checked 
against the Feynman rules. Two or three, exter- 
nal lines arc joined by means of a triple or quartic, 
vertex of the implemented model. In this way an 
off shell particle, which is represented by a spinor 
or polarization vector, is created. The off shell 
particles and/or external particles are joined in 
this way until the two parts of a considered topol- 
ogy are completely covered. If they match into a 
propagator then the topology is accepted. Once 
the topology has been accepted, the 'longer' part 
of it is further divided so that the Feynman dia- 
gram is made of 3 or 4 parts, joint to form a triple 
or quartic vertex of the model. Particles defined 
for one Feynman diagram are used as building 
blocks of other Feynman diagrams. The number 



of building blocks generated in carlomat is usu- 
ally smaller than in MadGraph. 

When the diagram is created, the correspond- 
ing 'particles' are used to construct the helicity 
amplitude, colour factor (matrix) and phase space 
parametrization which are stored on the disk. 
Once all the topologies have been checked subrou- 
tines for calculating the matrix element, colour 
matrix and phase space integration are written. 

2.3. Helicity amplitudes 

The helicity amplitudes are computed using the 
routines developed for MC programs ee4f gamma 
and eett6f [12 which have been improved 
and tailored to meet needs of the automatic gen- 
eration of amplitudes. In order to speed up the 
computation the MC summing over helicities has 
been implemented in the program. However, ex- 
plicit summing over helicities is also possible. 
While doing so, spinors or polarization vectors 
representing particles, both on- and off-shell ones, 
are computed only once, for all the helicities of the 
external particles they are made of, and stored in 
arrays, which is a novel feature with respect to 
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other programs, e.g. MadGraph. 

Possible poles in the propagators of unstable 
particles are regularized by constant widths which 
are introduced through the complex mass param- 
eters 

M| = m% - im B T B , B = W,Z, H, 
M t = yJm?-im t T t , (7) 

which replace masses in the corresponding prop- 
agators 

A ^ q) = q * - w v — ' 

A < \ 1 Q ( \ j + Mt 

Ap{q) = -?~W H SF{q) = ^Mf> 

both in the s- and i-channel Feynman diagrams. 
Propagators of a photon and gluon are taken in 
the Feynman gauge. The EW mixing parameter 
may be defined either real, referred to as the fixed 
width scheme (FWS), or complex, referred to as 
the complex mass scheme (CMS) 

sm Uw = i t~ or sm "w t- 

m% M| 

The colour matrix is calculated only once at the 
beginning of execution of the numerical program 
after having reduced its size with the use of the 
SU(3) algebra properties. 

2.4. Phase space integration 

A dedicated phase space parametrization is 
generated for each Feynman diagram. Mappings 
of the Breit-Wigner shape of the propagators of 
unstable particles and ~ 1/s behaviour of the 
photon and gluon propagators are performed. 
Dedicated treatment of soft and collinear exter- 
nal photons, as well as i-channel photon/gluon 
exchange is envisaged. 

The phase space parametrizations are incorpo- 
rated into a multichannel MC integration rou- 
tine for calculating total and differential cross sec- 
tions. The integration is performed iteratively. It 
starts with equal weights for all the kinematical 
channels and the weights with which each kine- 
matical channel contributes to the integral are 
determined anew after every iteration, carlomat 
can be used as MC generator of unweighted events 
as well. 



3. TESTS 

Matrix elements of many reactions with 6 par- 
ticles and several reactions with 7 particles in the 
final state have been checked against MadGraph 
for randomly selected sets of momenta. An agree- 
ment better than 13 digits has been found. 

Total cross sections of reactions e + e~ — > 4 
fermions and e + e~ — > 4 fermions and a photon 
have been checked against ee4f gamma [11] , and 
of e + e~ — * 6 fermions, relevant for the top quark 
pair production and decay, have been compared 
with eett6f . The results have agreed within one 
standard deviation of the MC integration. 

Moreover, checks against other MC programs 
have been made. Cross sections of the following 
top quark pair production reactions at the ILC: 

e + e~ — ► bbuddu, (8) 

e + e~ — » bbude~ v ei (9) 

e + e~ — ► bbv e e + fi~~ u^, (10) 

e + e~ — > bbv^n + n~ (11) 

computed with carlomat are compared against of 
the corresponding results of HELAC/PHEGAS and 
AMEGIC++ [T3] in Table [Q The cuts and initial 
parameters are those of [13] and, as in [13], the 
cross sections have been calculated with 10 6 calls 
in the MC integration, before the cuts have been 
applied. Satisfactory agreement also for all the 
other cross sections presented in [TJ] for about 
80 reactions has been found. However, for some 
reactions containing gluons, or e + e~ pair in the 
final state, more calls had to be used in carlomat 
in order to match the corresponding precision of 
[IB] , as the appropriate mappings have not yet 
been implemented in the program. Cross sections 
of reaction ([4]) without the gluon exchange con- 
tributions computed with carlomat and Whizard 
are compared in Table [2] Again a satisfactory 
agreement can be seen. 

4. SAMPLE RESULTS 

The capability of carlomat to handle multipar- 
ticle reactions with large numbers of the Feyn- 
man diagrams was demonstrated in [14j . where 
the off resonance background effects in the asso- 
ciated top quark pair and Higgs boson production 



Current status of carlomat, a program for automatic computation of lowest order cross sections 



5 



Table 1 

Cross sections in fb of ©-([Til) at y/s = 360 GeV 
(first row) and y/s = 500 GeV (second row). The 
cuts and initial parameters are as those of [13"] . 
The numbers in parentheses show the uncertainty 
of the last decimals. 



Reac. 


carlomat 


AMAGIC++ 


HELAC 


(8) 


32.98(11) 
50.31(19) 


32.90(15) 
49.74(21) 


33.05(14) 
50.20(13) 


© 


11.448(26) 
17.424(56) 


11.460(36) 
17.486(66) 


11.488(15) 
17.492(41) 




3.843(5) 
5.856(11) 


3.847(15) 
5.865(24) 


3.848(7) 
5.868(10) 


CD 


3.837(5) 
5.834(10) 


3.808(16) 
5.840(30) 


3.861(19) 
5.839(12) 



Table 2 

Cross sections in ab of (UJ) without the gluon ex- 
change contributions calculated with carlomat 
and Whizard. The cuts and initial parameters 
are as those of [14]. The numbers in parentheses 
show the uncertainty of the last decimal. 

y/s [GeV] carlomat Whizard 



500 
800 
1000 
2000 



7.80(3) 
66.8(1) 
61.4(1) 
28.5(1) 



7.76(2) 
67.3(1) 
61.8(2) 
28.1(3) 



at the linear e + e~ collider were studied. In [T4] . 
the effects were shown in reactions <j3j (El) - Here, 
we will address the question of the off resonance 
background contributions to the associated top 
quark pair and Higgs boson production in reac- 
tion © which is the hadronic detection channel 
of ([1]) that, in the unitary gauge and neglecting 
the Yukawa couplings to light fermions, receives 
contributions from 185 074 lowest order Feynman 
diagrams. 

As in [14 , let us identify jets with their orig- 
inal partons and define the following cuts on an 



angle between a quark and a beam and an angle 
between two quarks 

5° < %,beam) < 175°, %,?')> 10° (12) 

and a cut on a quark energy 

E q > 15 GeV, (13) 

which should allow to detect events with 8 sep- 
arate jets. Moreover, in order to reconstruct W 
bosons, top quarks and the Higgs boson let us as- 
sume 100% efficiency of b tagging and impose the 
following invariant mass cuts: 
a cut on the invariant mass of two non b jets 



60 GeV < 



(p~6 1 +P~b 2 Y 



1/2 



< 90 GeV, 



(14) 



a cut on the invariant mass of a & jet, b\, and two 
non b jets, b^, &~b 2 i 



(pb t +p~b 2 y 



1/2 



m t 



< 30 GeV (15) 



and an invariance mass cut on two b jets, &3 and 

& 4 , 



(Pb 



Pbi 



1/2 



< m 



bb > 



(16) 



with m c b f = 20 GeV, 5 GeV, 1 GeV. 

The total cross section of reaction (UJ) calcu- 
lated with the complete set of the lowest order 
Feynmamn diagrams, <7 a u, and with the 20 signal 
diagrams, representatives of which are depicted 
in Fig. [TJ Csigo with cuts (fl"2")) - (fl"6]) are shown in 
Table [3] The initial parameters used are those 
of [14] . We see that invariant mass cut (fT6]) . 
which has been imposed in order to identify the 
bb quark pair coming from the Higgs boson de- 
cay, reduces the off resonance background very 
efficiently while it practically does not alter the 
signal cross section. 

5. CONCLUSIONS AND OUTLOOK 

carlomat can be used for automatic computa- 
tion of cross sections of multiparticle reactions as 
it has been demonstrated in [14] and in Section 4. 
It can be used as an MC generator of unweighted 
events, too. 
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Table 3 

The cross sections a a \\ and cr s i g . of ([6]) for 3 differ- 
ent values of the invariant mass cut m^* of (|16|) . 
The other cuts are given by (fT^ - lfl!)]) and initial 
parameters used are as those of [2]. The num- 
bers in parentheses show the uncertainty of the 
last decimal. 



yfi [GeV] 


mtf [GeV] 


fail [ab] 


CT sig . [ab] 


500 


20 


13.46(5) 


8.72(1) 




5 


10.12(4) 


8.70(1) 




1 


8.98(4) 


8.67(1) 


800 


20 


164.5(4) 


128.6(1) 




5 


139.5(4) 


128.2(1) 




1 


129.6(2) 


127.7(1) 


1000 


20 


137.9(3) 


109.2(1) 




5 


117.7(5) 


109.1(1) 




1 


110.6(5) 


108.6(1) 


2000 


20 


44.2(2) 


36.17(4) 




5 


38.2(2) 


36.27(4) 




1 


36.5(1) 


36.14(4) 



In spite of all the successful checks presented 
in Section 3, further thorough comparisons with 
other existing MC generators still should be done. 

Interfaces to PDF's, or ISR within the struc- 
ture function approach are practically ready. In- 
terfaces to parton shower and hadronization pro- 
grams should be worked on. Extensions of SM 
can be implemented and the corresponding low- 
est order cross sections can be calculated in a fully 
automatic way. 

Leading SM radiative corrections can be im- 
plemented, if corresponding subroutines are pro- 
vided as it was done, e.g. in [15j . where factor- 
izablc EW 0(a) corrections were included for re- 
actions e + e~ — ► 6 fcrmions relevant for the top 
quark pair production and decay at a linear e + e~ 
collider. 
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